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Abstract 

Based on Novotny's transfer-matrix method, we simulated the (stacked) triangular Ising anti- 
ferromagnet embedded in the space with the dimensions variable in the range 2 < d < 3. Our aim 
is to investigate the criticality of the XY universality class for 2 < d < 3. For that purpose, we 
employed an extended version of the finite-size-scaling analysis developed by Novotny, who utilized 
this scheme to survey the Ising criticality (ferromagnet) for 1 < d < 3. Diagonalizing the transfer 
matrix for the system sizes N up to N = 17, we calculated the (i-dependent correlation-length 
critical exponent v (d). Our simulation result v(d) appears to interpolate smoothly the known two 
limiting cases, namely, the KT and d = 3 XY universality classes, and the intermediate behav- 
ior bears close resemblance to that of the analytical formula via the 1 /^-expansion technique. 
Methodological details including the modifications specific to the present model are reported. 



I. INTRODUCTION 



In analytical approaches, the spatial dimension d is treated as a continuously variable 
parameter, and correspondingly, various quantities such as the critical indexes are expressed 
explicitly in terms of the parameter d. Such an approach allows us to see how the criticality 
changes from the classical (meanfield like) one as the spatial dimension deviates from an 
either lower or upper critical dimension gradually. However, it is not quite obvious that such 
an analytical formula could be justified (realized) by actual first-principles simulations (and 
hopefully by experiments) with respect to realistic lattice models. In fact, in conventional 
computer-simulation approaches, one has to fix the (embedding) spatial dimension to a 
certain integral value, and thus the analysis on criticality has been restricted to the integral 
values of d inevitably. 

An attempt to circumvent such a restriction was made by Novotny [1-4]. His approach 
stems on a very formal expression for the transfer matrix so that the embedding spatial 
dimension can be varied continuously. (We explain his method in the next section. As 
anticipated naturally, this method is also of use in studying high dimensional (d > 3) 
systems. We refer readers to the references [1, 5] for this development.) Based on this 
formulation, he performed an extensive computer simulation, and surveyed the criticality of 
the Ising ferromagnet for 1 < d < 3. Astonishingly enough, he found that the numerical 
result is well described by both the 1 + e and 4 — e expansion formulas. In other words, 
his result clarifies that the analytical formulas for fractional values of d are meaningful 
in the sense that they are reproduced by the first-principles-simulation scheme. (Strictly 
speaking, he utilized two distinctive approaches to control the embedding spatial dimension. 
In Ref. [2], he varies the "connectivity" of the lattice, whereas in Refs. [3, 4], he twists the 
boundary condition to control the magnetic-do main- wall undulations.) Here, we stress that 
the Novotny approach is not a mere dimensional interpolation (crossover) that has been 
studied extensively in the past study [6, 7]. 

In this paper, we apply Novotny's method to the "stacked" triangular Ising antiferromag- 
net [8-21] embedded in the space with the variable dimensions 2 < d < 3. Because of its Z 6 
invariance, the model should exhibit the XY universality class at the magnetic transition 
point [22-26]. Our aim is to examine whether his method is applicable to generic problems 
other than the Ising universality class. We calculate the ^-dependent correlation-length 
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critical exponent. Thereby, we will show that the simulation result is comparable with the 
analytical 1/iV-expansion result up to 0(1/N) [27-30]. 

In fairness, it has to be mentioned that the critical phenomena for non-integer dimensions 
were studied extensively in the past [31-35]. In these works, the authors set up their lattice 
models on a fractal structure (the Sierpinski gasket) in order to realize a magnetism in the 
fractional spatial dimensions. Here, we stress that in our approach, the spatial dimension 
can be varied continuously within a range. 

The rest of this paper is organized as follows. In the next section, we explain how we 
constructed the transfer matrix for the stacked triangular antiferromagnet in 2 < d < 3. 
In Sec. Ill, we present the numerical results. Managing an extended finite-size scaling 
analysis, we estimate the correlation- length critical exponent in the range 2 < d < 3. In the 
last section, we present summary and discussions. 

II. CONSTRUCTION OF THE TRANSFER MATRIX FOR THE "STACKED" 
TRIANGULAR ANTIFERROMAGNET IN 2 < d < 3 

In this section, we set up the transfer matrix formalism to simulate the "stacked" tri- 
angular antiferromagnet embedded in the dimensions 2 < d < 3. Our formalism is based 
on Novotny's idea [1], with which he studied the Ising ferromagnet on (hyper) cubic lat- 
tices; in our preceding paper [5], we adopted his idea to the d = 3 Ising ferromagnet with 
plaquette-type interactions. Before going into details, we first set up a basis of our scheme; 
namely, we dwell on the particular d = 3 case. Then, we extend this preliminary basis to 
incorporate the embedding-spatial-dimension variation. We will also provide a number of 
technical modifications to improve the efficiency of the numerical simulation. 

We decompose the transfer matrix into the following two contributions, 

T = T i 0T l , (1) 

where the symbol denotes the Hadamard (element by element) matrix multiplication; 
note that the multiplication of local Boltzmann factors yields the global Boltzmann weight. 
As explained below (see also Fig. 1 for the geometrical structure of our finite-size cluster), 
the decomposed parts T 1 - and T" account for the contributions from the intra and inter 
plane (triangular lattice) interactions, respectively. 
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First, let us consider the component TK The matrix elements are given by the formula, 

T \\ _ U\A\A\ _ w S(j,l)S(j,2) w S(j,2)S(j,3) w S(j,N)SVA) (?] 
J-ij - \ l \^\J) - VV S(i,l)S(i,2) VV S(i,2)S(i,3) W S(i,N)S(i,l) ' \ Z ) 

where the indexes % and j specify the Ising spin configuration for both sides of the transfer- 
matrix slice; see Fig. 1. More specifically, we consider N spins for the transfer-matrix slice, 
and the index i denotes a spin configuration {S(i, 1), . . . , S(i, N)} arranged along the "leg." 
The factor W^g* stands for the local Boltzmann weight for a unit cell of the triangular 
lattice with the corner spins {Si, S 2 , S 3 , S 4 }. Explicitly, it is given by the following form, 

1 fJ 



S = exp 



r s 2 (S 1 S 2 + S 2 S 4 + S 4 S 3 + S 3 S ± ) + JS 1 S i 



(3) 



(The denominator of the coupling constant is intended to avoid double counting.) Here, 
the parameter T denotes the temperature, and the parameter J stands for the intra-plane 
antiferromagnetic interaction constant. Hereafter, we choose J as the unit of energy; namely, 
we set J = 1. It is to be noted that the component TH (with T L ignored) leads to the transfer- 
matrix for a sheet of triangular antiferromagnet. In other words, the remaining component 
T 1 - should raise the dimensionality to d — 3 through introducing the inter-plane interactions. 
This is an essential idea of the Novotny method [2]. 

Second, we consider the component T 1 - that accounts for the inter-plane interaction. The 
explicit matrix elements are given by the following formula, 

T£ = (i\BP v \i), (4) 

with the interaction distance v. The matrix B is given by the formula, 

= ^s(i,i)s(j,i)^s(i,2)s(j,2) ' ' ' Ws(i,3)S0'>3)' (5) 

with Wg lS2 = exp(—jSiS 2 /T) and the inter-plane interaction j. The matrix P denotes 
the translational operator: That is, with one operation of P, a spin arrangement {S(i,m)} 
shifts to {S(i,m + 1)}; note that the periodic boundary condition is imposed. An explicit 
representation of P is given in our preceding paper [5]. Because of the insertion of P v , 
the interaction B bridges the i>th neighbor pairs along the leg, and so it brings about the 
desired inter-plane interactions. As a matter of fact, in Fig. 1, we notice that the alignment 
of spins is folded into a rectangular shape with the edge lengths v and N/v. It is an essential 
idea of Novotny that the operation P v is still meaningful, although the power v is not an 



integral value. This rather remarkable fact renders freedom that one can construct the d = 3 
transfer-matrix systematically with arbitrary number of spins N. 

Based on the above formalism, we readily simulate the stacked triangular antiferromagnet 
in d = 3. In the following, we propose a scheme to tune the embedding spatial dimension 
continuously. Moreover, we will also provide a number of technical modifications, aiming to 
improve the efficiency of the simulation. 

There are two controllable parameters for the dimension variation. That is, the inter- 
plane interaction j and the interaction distance v: Apparently, the limit j — > reduces the 
system to a sheet of triangular lattice. On the other hand, for large v, the stack width N/v 
decreases, and eventually, at v — N, the system reduces to a sheet of triangular lattice as 
well; note the identity P N = 1 owing to the periodic boundary condition. In this paper, 
we adopt the former scheme. Namely, we will tune the parameter j, fixing the interaction 
distance v to a moderate value v = 0.27N. This choice is based on our observation that the 
finite-size-scaling behaviors become quite systematic for v ~ N/n (n: integer), particularly 
at n — 4. 

Lastly, let us explain a number of technical modifications to improve the efficiency of the 
simulation: We propose the following replacement, 

T± = (i\BP»\i) -> T$(v) = (i\Bpv\i)(i\p-vB\i). (6) 

(Note that correspondingly, we need to replace the temperature T with 2T in order to 
compensate the duplication.) With this trick, the transfer-matrix elements become real: 
Otherwise, the elements are complex for even values of N. As a matter of fact, in the past 
simulations [1, 2, 5], those cases of even values of N were excluded. Such an exclusion is 
obviously disadvantageous in the subsequent data analysis, because the available systems 
sizes are restricted severely. In our simulation, because of the above trick, we are able to 
consider arbitrary system sizes. In addition to this, we symmetrize the transfer matrix [2] 
with the following replacement, 

T^^T^QT^-v). (7) 

(Similarly as the above, we need to re-define the temperature T — > 2T.) With this sym- 
metrization, the symmetry of the descending (j — N, N — 1, . . . ) and ascending (j = 1, 2, . . . ) 
directions along the leg become completely restored. We observed a significant improvement 
of the finite-size scaling behavior due to this symmetrization. 
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III. NUMERICAL RESULTS 



In the preceding section, we developed a transfer-matrix formalism for the stacked tri- 
angular antiferromagnet embedded in the fractional spatial dimensions 2 < d < 3. In this 
section, we present the numerical results calculated by means of the exact-diagonalization 
method for the system sizes up to N = 17. We analyze the data with the extended finite-size 
scaling analysis [2] which allows us to estimate the "effective" dimension <i e //- The effective 
dimension plays a significant role in the subsequent analysis of the criticality of the magnetic 
transition. 

A. Effective dimension: Extended finite-size scaling analysis [2] 

Before going into detailed analysis on the criticality, we need to estimate the effective 
dimension d e ff [2]: At the critical point, the correlation length should be comparable to 
the linear dimension of the finite cluster. Hence, the correlation length £ should obey the 
formula £ ~ ^/(^ff- 1 )^ no te that a transfer-matrix slice contains N lattice points, and 
its embedding spatial dimension should be <i e // — 1- This formula immediately yields an 
estimate for the effective dimension, 



for a pair of system sizes (N, N') . By means of the transfer- matrix method, the correlation 
length is calculated immediately: Using the largest and the next-largest eigenvalues, namely, 
Ai and A 2 of the transfer matrix, we obtain the correlation length as £ = l/ln(Ai/A 2 ). 
Provided by this, according to Ref. [2], we are able to determine both the critical temperature 
T c (Ni, N 2 ; N 3 , N4) and the effective dimension d e ff(Ni, N 2 ; N 3 , JV4) so that they satisfy the 
following equation, 



d eff (N u N 2 ; N 3 , N 4 ) = d"y a (T c (N u N 2 ; N 3 , A 4 )) = d^f^T^ N 2 ; N 3 , N 4 )), (9) 



for the set of the system sizes (Ni, N 2 ; N 3 , N4). To summarize, in the extended finite- 
size scaling analysis, the spatial dimension is not a given constant, but a parameter that 
is to be determined a posteriori with the data analysis of the correlation length £. As 
mentioned in the above, in the transfer-matrix method, the correlation length is calculated 
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quite straightforwardly. In that sense, the transfer-matrix approach is suitable to this type 
of finite-size-scaling analysis. 

In order to examine the validity of the scaling parameters, d e ff and T c , determined with 
the above method, we plotted, in Fig. 2, the scaled correlation length (T — T^)N 1 ^ v ^ de ^~ 1 ^ > - 
f/jVVte/z-i) for d eff = 2.50, T c = 5.13, l/v = 1.36, j = 1.5, and TV = 13, 14, . . . , 17. 
These scaling parameters, namely, d e ff and T c , are determined from the set of system sizes 
(14, 16; 13, 16) via the extended finite-size-scaling analysis. (We explain how we determined 
l/v — 1.36 afterward.) We see that the scaled data collapse into a scaling function quite 
satisfactorily. Hence, we confirm that the scaling parameters, d e // = 2.50 and T c = 5.13, 
are indeed meaningful. More significantly, we stress that our simulation data should be 
described under the assumption that the effective dimension takes such a fractional value. 

In order to analyze the criticality further in detail, we calculated the Roomany-Wyld 
approximative (3 function, which is given by the following formula [36], 

In Fig. 3, we plotted (3\±,iq(T) for the same parameters as those of Fig. 2. The functional 
form of this (3 function seems to be almost straight, indicating that the corrections to the 
finite-size scaling are almost negligible. This fact accounts for the good data collapse of Fig. 
2 shown above. From the slope of this (3 function at the transition point T = T c , we are 
able to estimate the inverse of the correlation-length critical exponent as 1/Vi4,i6 = 1-36. To 
summarize, we estimated the exponent 1/Vjv,jv' from the set of system sizes (TV, TV'). In prior 
to this analysis, we should determine the scaling parameters, d e ff and T c , from the extended 
finite size scaling analysis for (TVi, TV 2 ; TV 3 , TV 4 ). We will exploit the <i e //-dependence of l/v 
in the next subsection. 

Lastly, we exploit the region in close vicinity of the lower critical dimension d — 2, at 
which the KT-type singularity should emerge. In Fig. 4, we present the function (3\a,ig{T) 
for j = 0.9, T c = 2.74, and d e ff = 2.07; these scaling parameters were determined from the 
set of system sizes (14, 16; 15, 17). In contrast to the behavior shown in Fig. 3, the f3 function 
is curved particularly in the vicinity of the transition point. Actually, we estimate the slope 
(critical exponent) as 1/Vi4,i6 = 0.804 which is considerably suppressed, compared with that 
of Fig. 3. This feature indicates that an essentially singular-type critical behavior emerges 
as d e ff — > 2. For T < T c , the (3 function starts to increase, and eventually, it becomes 



even positive in the low-temperature regime. Such a feature may reflect an instability to 
the Z 6 -symmetry-broken phase. Actually, it has been known [22] that right at d — 2, an 
additional phase transition of the KT type takes place at a low temperature, where the 
Z 6 -symmetry-breaking field becomes marginally relevant. The simulation data around this 
regime may be affected by the notorious logarithmic corrections to the finite-size-scaling 
behavior, that are inherent to the KT-type critical phenomenon. 

B. Correlation-length critical exponent l/z/(d e //) for d e ff 

In the above, we analyzed the criticality of the magnetic transition for j = 1.5 and 0.9 
in terms of the effective dimension d e ff. Managing the similar analysis for various j, we are 
able to survey the d e ff- dependence of the critical exponent l/v(d e ff ). In Fig. 5, we plotted 
the critical exponent 1/V for various d e ff. The exponent was determined from the set of 
system sizes, (■) (14, 16), (x) (15, 17), (*) (14, 16), and (□) (13, 17), for respective symbols, 
and the corresponding scaling parameters, namely, d e ff and T c , had been determined from 
the set of system sizes (■) (14, 16; 15, 17), (x) (15, 17; 13, 17), (*) (14, 16; 13, 16), and (□) 
(15, 17; 13, 17), respectively. We also plotted a result \jv = 1.48909(60) for the d = 3 XY 
universality class [37] with the symbol +. We notice that our numerical results interpolate 
smoothly the known limiting cases of KT (1/V = 0) and d = 3 XY universality classes. As 
for a comparison, with a dotted line, we presented the 1/iV-expansion-approximation result 



with N = 2 and e = d — 2. Rather remarkably, we see that the simulation data and 
the 1/iV-expansion result exhibit similar intermediate behaviors. In other words, we can 
make contact with such a dimensional-regularized analytical expression via the computer 
simulation calculation. On closer inspection, however, it seems that our first-principles 
simulation predicts even more convex-like functional form. Actually, our simulation result 
suggests a notable steep increase around d = 2. 

It is to be noted that our data extend to the regime exceeding the threshold d = 3. As 
a matter of fact, we intended to cover the parameter range 2 < d < 3, when we constructed 
the transfer matrix in Sec. II. However, such a feature was also observed in the preceding 



up to 0(1/N) [27-30], 
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study of the Ising ferromagnet [3]. In the study, the author reported that the effective 
dimension does exceed the intended range, and astonishingly enough, the data are still in 
good agreement with the (4 — e)-expansion result. We expect that our result makes sense 
even for d e ff > 3 as well. Our data seem to approach to the meanfield value \jv — 2 as 
d — > 4 rather directly than the 1/iV-expansion- approximation result. 

In close vicinity of the lower critical dimension, particularly, around d e ff ~ 2.3, the 
numerical data turn out to be scattered. As noted in the preceding subsection, this data 
scatter should be attributed to the logarithmic corrections to the finite-size scaling, which 
are inherent to the KT-type critical behavior. 

IV. SUMMARY AND DISCUSSIONS 

In this paper, we developed a transfer-matrix scheme to simulate the "stacked" triangular 
antiferromagnet embedded in 2 < d < 3. Our scheme is based on Novotny's idea, which has 
been applied to the Ising ferromagnet (universality class) in 1 < d < 3 successfully [2-4]. 
Here, we studied the XY universality class in 2 < d < 3; the triangular antiferromagnet 
should belong to the XY universality class due to the Z e symmetry [22-26] . The numerical 
data are analyzed in terms of the extended finite-size-scaling method [2] , which allows us to 
estimate the effective dimension d e ff. Thereby, we obtained the o? e //-dependent correlation- 
length critical exponent l/v(d e ff); see Fig. 5. We notice that our first-principles data 
interpolate smoothly the known limiting cases of both KT and d = 3 XY universality 
classes. Furthermore, we found that the intermediate behavior bears close resemblance 
to that of the analytical formula via the 1/A-expansion technique. In other wards, by 
means of the computer simulation method, we are able to check (support) the validity of 
the dimensional- regularized formulas with the e and 1/N expansion techniques. On closer 
inspection, our simulation result suggests a notable steep increase around d = 2. 

The present study on the stacked triangular antiferromagnet shows that the Novotny 
method would be generic, and it should be applicable to a wide variety of universality classes 
other than the Ising universality. So far, the transfer-matrix approach has been restricted 
to the problems in two dimensions, because it requires huge computer memory space as the 
system size increases. (Although the density-matrix-renormalization-group method resolves 
this difficulty to a considerable extent, its extension to d — 3 is still a current topics underway 
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[38].) With the aid of the Novotny method, we are able to construct the transfer matrix 
for d > 2 quite systematically with modest (actually, arbitrary) number of constituent 
spins. Moreover, we can survey the criticality even in the fractional dimensions by means of 
the extended finite-size-scaling analysis. This opens a way to re-examine the longstanding 
problems in three dimensions such as the chiral universality [39] and the Lifshitz-type multi- 
critical phenomenon [40]. 

One may wonder what affects the variation of d e ff most significantly: Actually, there 
have been known two approaches in order to control d e ff in the past studies of the Ising 
ferromagnet. In Ref. [2], the interaction distance v, in other words, the connectivity of 
the finite-size cluster, is tuned carefully. On the other hand, in Refs. [3, 4], the boundary 
condition is twisted so as to control the thermal undulations of the magnetic domain walls. 
In our scheme, as explained in Sec. II, we fixed the lattice connectivity v, and rather 
varied the inter-plane interaction j: In this sense, we took an advantage that our system 
(stacked triangular antiferromagnet) is, by nature, spatially anisotropic, and so we are able 
to tune the inter-plane interaction freely. We suspect that our case may belong to the 
latter category. That is, the inter-plane interaction, somehow, controls the domain-wall 
undulations effectively. This interpretation is based on the fact that our system is subjected 
to the magnetic frustration due to the triangular antiferromagnetism, and the domain walls 
should be created inevitably. It is desirable that this mechanism would be exploited in the 
future study. 
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FIG. 2: Scaling plot of the correlation length (T-T c )ArVKd e //-i)_f/jvWe//-i) j s s hown for j = 1.5 
and the system sizes AT = 13, 14, ... , 17. The symbols +, x, *, □, and ■ denote the system sizes 
of N = 13, 14, 15, 16, and 17, respectively. The scaling parameters are set to be 1/v = 1.36, 
d e ff = 2.50, and T c = 5.13; see text for details. We see that a good data collapse is achieved under 
a fractional value of the effective dimension d e // = 2.50. 




FIG. 3: We plotted the beta function, /3u : ie{T) (10), for the same parameters as Fig. 2. From the 
slope at the transition point, we estimate the correlationdength critical exponent as 1/v = 1.36. 
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FIG. 4: The beta function, /3 Ujl6 (T) (10), is plotted for j = 0.9 and d eff = 2.07. We notice that 
for d e jf « 2, the beta function gets curved, indicating that a nonstandard (essentially singular) 
criticality emerges. From the slope at the transition point, we estimate the correlation-length 
critical exponent as 1/u = 0.804. 
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FIG. 5: The inverse of the correlation-length critical exponent l/u(d e ff) is plotted for the effective 
dimension d e //. The exponent is determined from the pair of system sizes, (■) (14, 16) under the 
extended scaling analysis [2] for (14, 16; 15, 17), (x) (15,17) with (15, 17; 13, 17), (*) (14,16) with 
(14, 16; 13, 16), and (□) (13,17) with (15, 17; 13, 17), for respective symbols; see text for details. 
We also plotted the 1/iV-expansion-approximation result up to 0(l/iV) [27-30]. The symbol + 
denotes a result for the d = 3 XY model [37]. 




14 



